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. We consider wormhole geometries subject to a gravitational action consisting of 

non-linear powers of the Ricci scalar. Specifically, wormhole throats are studied 
in the case where Einstein gravity is supplemented with a Ricci-squared and 
inverse Ricci term. In this modified theory it is found that static wormhole 
q-( throats respecting the weak energy condition can exist. The analysis is done 

locally in the vicinity of the throat, which eliminates certain restrictions on the 
models introduced by considering the global topology. 
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1 Introduction 

Wormhole type solutions in Einstein gravity have long been studied. Flamm pQ seems 
to have been the first to study such objects while considering the then newly discov- 
ered solution of Schwarzschild. In 1935, Einstein and Rosen considered wormhole type 
bridges as potential models for elementary particles |2j: the famous Einstein- Rosen 
bridge. The field lay dormant for approximately two decades until the consideration by 
Wheeler of the possibility of a space-time foam. This foam, due to violent fluctuations 
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in the metric when considering quantum gravity, could be viewed as wormhole-like 
structures permeating through space-time which would potentially be manifest at en- 
ergies near the Planck scale |3,. As well, wormholes have relevance to issues such as 
chronology protection 0] , topology change |5J and in studies of horizons and singular- 
ities 0, [7j. 

There has been a sizeable amount of literature produced over the past 15 years 
on the issue of Lorentzian as well as Euclidean wormholes. This re-newed interest in 
the subject was mainly sparked by the work of Morris and Thorne |S] where they 
considered static, traversable wormholes for pedagogical purposes. In their work, and 
subsequent papers by others, it was found that a static wormhole throat could not 
be supported unless the weak energy condition (WEC) was violated. The violation 
could be confined to a small region in the vicinity of the throat but it must occur to 
prevent the throat from collapsing (see [§] , [ID] , [U] , [J2] , |Hj and references therein) . 
Non-spherical geometries have also been studied (some examples are [H], (THj, [To] . 
|17j). For a more complete list of pre-1996 references, as well as an excellent exposition 
on the subject, the reader is referred to the book by Visser [T5j . 

There have been many studies of static wormhole solutions employing various 
exotic matter models within Einstein gravity. However, it is possible that the action 
governing gravitational dynamics is not the Einstein-Hilbert one, but consists of a 
more complicated Lagrangian, which reduces to the usual Einstein-Hilbert action in 
some limit. In this vein theories with generalized actions of the form : 



with C(R) being some function of the Ricci scalar, have been employed to explain 
various phenomena [20] . The most popular of these supplements the Einstein-Hilbert 
Lagrangian with an R 2 term so that £(R) = R + aR 2 (for recent work the reader is 
referred to |2*T| ) . This leads to modifications which could be utilized to drive inflation 
purely from the gravitational sector (the Starobinsky inflationary theory 22j). A 
thorough study of solutions in R m cosmology may be found in the recent papers [22] 
as well as [21]. 

A more recent modification which has received attention is the addition of a 
1/R term in the Lagrangian. The major motivation for this modification is in the 
realm of matter dominated era cosmology. For example, it is thought that perhaps 
such a term may contribute to the present day acceleration of the universe without the 
need for an exotic dark energy J2H], |26j . |27j . [2*5] , [2*U] . Inverse Ricci terms may also 
appear in certain sectors of string/M theory (see [HO], [HI] for example). These are 
interesting ideas although it has been shown that adding inverse powers of the Ricci 
scalar to the gravitational Lagrangian is accompanied by some potential problems [32*] . 




(1) 



m, a m- 



Notations and conventions here follow those of 
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Several issues are those of stability and the Newtonian limit. It has been 
argued that, in both the metric variation and the Palatini formulation of the 1/R 
theory, instabilities or potentially unphysical weak- field limits may exist. The Newto- 
nian limit for the metric variation theory has been carefully studied by Dick [37] who 
derived criteria for well defined Newtonian limit. (A similar analysis has been done by 
Dominguez and Barraco [33] in the Palatini formulation). We address this, and the 
issue of stability, later in the paper. It has also been shown by Chiba [33] that if the 
1/R sector is to drive the present day cosmological acceleration, the theory will yield 
results incompatible with solar system experiments. However, this does not preclude 
singular non-linear gravity theories which do not dominate present day cosmological 
evolution yet whose effects may be important under conditions as those potentially 
found near wormhole throats. This is also the motivation in a recent paper [33] where 
the future evolution of the universe is considered with inverse curvature Lagrangians. 

For the above reasons we believe it is of interest to study wormhole geometries 
in a theory whose action is of the form (|1J). That is, can the modified gravitational 
sector of the theory eliminate the need for exotic matter in supporting a static worm- 
hole throat? The most complete case for which the analysis is tractable seems to be 
the choice C(R) = C-iR^ 1 + R + C2R 2 , the c n 's being the "coupling" constants of the 
theory. The R 2 theory was considered by Ghoroku and Soma |3E] • in their meticulous 
study it was concluded that, under the assumption that an asymptotically flat global 
solution exists, no WEC respecting wormhole can exist in such a theory. The 1/R 
theory has not been considered in the context of wormholes. 

Here we study wormholes from a local geometric perspective, without con- 
straining the asymptotics. As pointed out in [3JJ], geometric information is less lim- 
iting than topology to the issue of defining and locating wormhole throats. We find 
that in the local analysis, R 2 terms permit the existence of a WEC respecting throat. 
The near throat solution can be patched, via appropriate junction conditions, to other 
solutions with desired far throat geometry and topology. This approach is often taken 
in, for example, studies of stellar structure or gravitational collapse. We briefly discuss 
the junction conditions in a later section. 

There are compelling reasons to believe that if gravitational fields are dictated 
by an action such as (0), then the non- linear contributions could be important in the 
case of wormholes. The R 2 terms would certainly be important in the high curvature 
regime which may be present in the vicinity of wormhole throats. As well, higher 
derivative contributions, introduced by both the R 2 and 1/R terms, may be important 
since the wormhole structure only places restrictions on the spatial components of the 
metric and its first derivatives. We show below that, in the WEC, the dominant terms 
near the throat are in fact those introduced by these non-linear terms, regardless of 
how small the coupling constants may be. It is found that the near throat region can 
respect the WEC in the higher derivative theory. 
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2 R m gravity and the wormhole geometry 

The general R m action is given by: 

S = ^[d i x^Yl c ™ Rm + S ™*"-> ( 2 ) 

m 

which, in the metric variation theory, gives rise to the equation of motion: 
W v mc m R m - x - Ut £ c m R m + SI J2 mc ™ ( Rm ^) ■* ~ E mCm ( i?m_1 ) ' 



m 



=:G*+Hl = teT* , (3a) 
with 

G%:=R%- X -RS% . (3b) 

In the subsequent text, we consider only m = — 1, 1 and 2. In such a theory, the vacua 
(T% = 0) are the de Sitter and anti-de Sitter solutions with effective cosmological 
constant, A, given by: 

A = ±^V^=f forc -K°- (4) 

Here, Rq is the Ricci scalar of the vacuum solution. 

In the curvature coordinates often employed in wormhole studies, the static 
line element may be written as: 

ds 2 = - e 7(r) dt 2 + e a{r) dr 2 + r 2 d6 2 + r 2 sin 2 d<f) 2 , (5) 
-oo <t<oo, < r < r < n, < 9 < ir, < < 2n . 

As we are interested in the region near the throat, we consider here only the local 
geometry in this vicinity. The t =constant, 9 = tt/2 profile curve is displayed in 
figure H where it is embedded in a higher dimensional space, whose extra coordinate 
is denoted x. Geometric structure is given to the wormhole manifold by imposing 
the constraint x = P±(r), with P±(r) being the shape function of the profile curve. 
(The + corresponds to the upper portion of the curve while the — corresponds to the 
lower portion). The wormhole is obtained by considering the surface of revolution 
generated when rotating the profile curve about the x-axis (inset). The lower half of 
the wormhole is obtained via a similar construction with a potentially different profile 
function. The line element (JSJ may be written is terms of the shape function as: 

d s 2 = - e J(r)± d f + {! + [P ± r ( r )] 2 } dr 2 + r 2 d9 2 + r 2 sin 2 9 d(j) 2 . (6) 

It is sufficient to study only the upper portion of the wormhole (P+(r)) and we therefore 
suppress the subscript in the following. 

The restrictions on P(r) are evident from the figure: 
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Figure 1: Profile curve, P±(r), in the t =constant and 6 = tt/2 submanifold. The wormhole 
is generated via rotation about the x-axis (inset). 

• P(r) must possess positive first derivative, at least in the vicinity of the throat. 

• P{r) must possess negative second derivative near the vicinity of the throat. 

• P, r {f) — ► +00 a s r — > r o (the throat). 

It is this last condition which is problematic in the study of wormhole throats. Al- 
though the infinities should cancel out in calculations of physical quantities involving 
non-singular throats, the analysis (both numerical and analytic) is hampered by the 
presence of these infinities. A resolution is to switch to a different coordinate system 
where such infinities are not present. 

Another coordinate system, sometimes utilized in wormhole analysis, is the 
the proper length gauge where the line element takes on the form: 

ds 2 = -e M dt 2 + dp 2 + r 2 (p) d6 2 + r 2 (p) sin 2 9 d<f . (7) 

The coordinate p represents the (signed) proper length coordinate in the radial direc- 
tion: 

p(r) = ± ^ e a±(r ' )/2 dr' . (8) 

Jr 

The throat in this chart is located at p — 0. Although continuous at the throat, the 
equations of motion are significantly more complicated in this gauge. There is also 
the problem that the profile function is no longer explicit in the metric which makes 
the mathematical wormhole engineering a more difficult task. 

In order to remedy these issues, we instead use a set of coordinates utilized 
in The construction is briefly summarized here: One begins with the standard 
curvature coordinates giving rise to figured! The profile function, x = P(r) is inverted 
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and instead we consider the new function given by r = Q{x) = P^ 1 {x). Essentially, 
we take the wormhole and rotate it as illustrated in figure El In this new gauge, the 
line element takes the form: 

ds 2 = -e x(x) dt 2 + {1 + [Q, x {x)] 2 } dx 2 + Q 2 {x) d6 2 + Q 2 (x) sin 2 6 d<f) 2 . (9) 

Notice that this has both the advantage of requiring a single coordinate chart as well 
as retaining the profile function in the metric. The infinite derivative at the throat is 
replaced with a vanishing derivative and therefore the metric singularity is eliminated, 
causing no trouble for analysis. Further, given the complexity of the equations of 
motion, we consider here the zero-tidal force class of wormholes as studied in [S], |^Tj . 
[12], and others. This amounts to setting X(x) constant, at least near the vicinity 
of the throat (actually, for the local analysis presented here, X(x) need not be constant 
so long as the first few derivatives of X(x) vanish at the wormhole throat as not to 
affect the analytic expansions presented below). 




Figure 2: Wormhole profile curve in the t =constant and 9 = 7r/2 sub manifold using the 
rotated system. The profile function is given by r = Q(x) = P~ 1 (x) and the radius of the 
throat is Qq. As before, the wormhole is generated via rotation about the x-axis (inset). 

Here the restrictions for a throat are: 

• The derivative of Q(x) must change sign at the throat. Q )X (x) < for x < 
and Q >x (x) > for x > 0. 

• Q(x) must possess positive second derivative near the vicinity of the throat. 

• Q,x( X )\x=0 = o. 

Save for the assumption of static spherical symmetry, the only other assump- 
tion we make regarding the spatial metric is that it is analytic in the throat region. A 
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general form of profile curve, suitable for studying any near throat geometry, is given 
by HO]: 

Q(x) =Q + A 2 x 2n e h{x \ (10a) 
Q tX {x) =A 2 x 2n e h(x) pnaT 1 + h x (x)] , (10b) 

Q, XjX =A 2 x 2n e h(x) { [2nx~ l + h x (x)] 2 + h, x , x (x) - 2m- 2 } . (10c) 

Here, A and Qo > are constants, n is a positive integer sufficiently large to pre- 
clude singularities and h(x) is an arbitrary analytic function of x. Higher derivatives 
are required for the field equations (}3"a|) but no restriction is placed on them by the 
requirement that the metric describe a wormhole throat. The presence of arbitrary 
functions allows the modeling of infintely many near-throat geometries. It is worth 
noting at this point that the above restrictions may be relaxed. One could, for exam- 
ple, leave Q(x) as completely arbitrary, save for the analyticity requirement and the 
assumptions on the signs of the first two derivatives. This yields recondite analysis 
but does not change the qualitative results of this paper. 



3 Weak energy condition 

The weak energy condition is the statement that no time-like observer measures a 
local energy density which is negative. It is usually formulated as: 

T^u" > V time-like u a , (11) 

the u a being the observer's four- velocity. The solution for the general WEC in spherical 
symmetry has been solved in [33] . For the static model presented here, it is sufficient 
to study the limiting trajectories in each of the three principal spatial directions. If the 
WEC is satisfied for these trajectories, it will be satisfied for all trajectories. Enforcing 
u a u a = — 1, the limiting cases of (|TT|) reduce to the four conditions: 

—T° >0 , (12a) 
— Tq + T\ >0, no summation over i . (12b) 

Using the metric ansatz (|10a|) in (|3a|) . we expand the components of the stress-energy 
tensor about the throat (x = 0). The results are summarized here: 

- T °o ^17)4 m + (1 + 96Q^VW) {Qlc^ + 8c*)] 

+ -^120A 2 e h{0 \h x ) x=0 (Q^c_! + 8c 2 ) x + ... for n = 2, (13a) 



„2n-4 



^( i+ ^ + 3) + ^ vB (*- + a 

+ . . . for n > 2, (13b) 
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along with: 



io + il ^(2n-4)! Ae [ QoC - 1+ Qo) 



2n-2 



X 



^(^^ (0, (M^o(Q^ + |)^ + ..., (14a) 
+ ... ^-r° + r 3 3 . (i4b) 

It is easy to see from the above expressions that the WEC, in principle, may be 
satisfied in the throat vicinity. The WEC expression most sensitive to deviations from 
Einstein gravity is (jl4aj) (this is the term whose near throat region must violate the 
WEC in Einstein gravity). 

The condition for well defined Newtonian limit in singular non-linear gravity 
models provided by Dick [37] read that, given i? > 0: 

\C(R)C"(R)\ {R=Ro ^l. (15) 

Here primes denote differentiation with respect to R and Ro is the Ricci scalar of the 
supported vacuum solution. It was shown that these conditions yield an acceptable 
Newtonain limit on length scales <C (i?o) -1 ^ 2 - The vacuum of (|5ajl possesses Ricci 
scalar given by (@J) (for m = —1, 1, 2 terms). Therefore, sufficient conditions for (|3a|) 
to possess correct Newtonian limit are: 

c_!<0, R = + y /3\c- 1 \, c 2 sa \ =. (16) 

3 J/ Vl c -il 

The first of these conditions guarantees the existence of a maximally symmetric grav- 
itational vacuum solution (de Sitter) whereas the second and third conditions satisfy 
Dick's restrictions [2Z1- Therefore, if (J16)) holds, WEC violation will not occur as long 
as the condition 

> Qo (17) 



V27|c_i| 3 / 2 

is satisfied with n > 1. For a theory with only Ricci squared terms the sole restriction 
is C2 > whereas near throat WEC violation is more severe in the case where only 
1/R modifications are present. 

Regarding the issue of stability, it is known that the theory governed by the 
Lagrangian in (0) with m — — 1, 1 and 2 is equivalent to, in the Einstein frame, 
gravitation coupled to a scalar field. The stability of such a theory has been studied 
in where effective potential techniques were utilized to analyze the stability. The 
effective potential, V(R), is given by 



(21 




| + c 2 R 3 ) c 2 R 3 


(R 2 + \ 




+ 2c 2 i? 3 ) 2 
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If the last term in (j!6|) is treated as an equality, the potential has real stationary 
points at R = +l/(3c 2 ), R = and R = — l/(3c 2 ) with the product 

\c-i\c\ = ^. (19) 

Under condition (|19|) the critical point at R = +l/(3c2) is a rising point of inflection 
whereas the critical point at R = — l/(3c2) is a minimum. However, the last criterion 
of (fT^j) is an approximate equality and it is easy to see that even a slight perturbation 
in the value of ()19|) will shift the inflection point and turn it into a minimum. Thus 
yielding both weak field stability and consistent weak field limit. We plot this region 
of the potential in figure El It should be noted that in strong gravitational fields (such 
as those potentially present near the wormhole throat), R need not be restricted to an 
interval near this value. However, this analysis is useful in illustrating that the theory 
possesses acceptable weak field behaviour. 

Effective Potential 

0.085485 H 



C_2V 

0.0854845 - 




0.34 0.35 
e_2 R 



Figure 3: Effective potential for stability analysis with |c_i |c| ~ 0.038. 



4 Junction conditions 

We briefly comment here on the patching at the junction away from the throat (located 
at x = > 0). As mentioned above, a patching of the near throat solution to a WEC 
respecting space-time could yield a WEC respecting wormhole. There are several 
physically acceptable junction conditions one could employ at the boundary (see, for 
example, [12], [HI], jH>], $T\, jlH] for applications to wormholes in Einstein gravity). 
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One condition is that of Synge [49J which reads: 

^ [T^% =Xt - e = ^T^h\ =Xt+e . (20) 

Here h u is the outward pointing unit normal vector of the junction surface. 

In Einstein gravity, another condition is the Israel-Sen-Lanczos-Darmois (ISLD) 
junction conditions j5U], [HI], [52] - [53] which reads: 

lim [K^x* - e)] = lim [Ku V {x* + e)] , (21) 

with K^y being the extrinsic curvature tensor of the junction hyper-surface. It has 
been shown that, on a time-like boundary in a spherically symmetric i?-domain, the 
above two conditions are equivalent in Einstein gravity jH]. 

The Synge junction condition (J20|) carries over, in a straight forward manner, 
into higher derivative gravity. The ISLD condition (|2*T|) . relying on first derivatives of 
the metric, needs to be modified to be valid in higher derivative theories. One way this 
could be accomplished is to consider the boundary terms which arise when applying 
the variational principle to The continuity of such terms at the junctions will 
lead to an appropriate junction condition. Without detailed analysis, it is probably 
safe to state that continuity of up to fourth derivatives of the metric will satisfy such 
a junction condition (although this restriction can most likely be relaxed somewhat). 
This guarantees continuity of G^ + H^ (or, equivalently here, T^). The above solution 
may be patched, for example, to an anisotropic fluid which is characterized by: 

T$ = (p + p±) u^u u + P± 5t + ( Pll - P± ) s<+s u , (22) 

with the fluid four-velocity, which is perpendicular to the vector s M . The quantities 
p, p^, pi! represent the energy density, transverse pressure and parallel pressure of the 
fluid respectively. Other boundaries further away from the throat, such as a matter- 
vacuum boundary, may similarly be considered. 

5 Concluding remarks 

In summary, wormhole throats were studied in a gravitational theory governed by the 
Einstein-Hilbert Lagrangian supplemented with 1/R and R 2 Ricci scalar terms. The 
resulting equations of motion was utilized to study energy conditions in the vicinity 
of the wormhole throat. It was found that the weak energy condition (WEC) may 
be respected in the throat region in the modified theory. The conditions for a WEC 
respecting throat are compatible with those required for stability and an acceptable 
Newtonian limit. Away from the throat, the system cold be joined to WEC respecting 
matter solutions and, at other junctions, even the vacuum. It would be of interest to 
study these junction conditions in detail in future work. Other curvature invariant 
contributions could also be investigated. 



10 



References 

L. Flamm, Phys. Z. 17 (1916) 448. 
A. Einstein and N. Rosen, Ann. Phys. 2 (1935) 242. 
J. A. Wheeler, Phys. Rev. 48 (1957) 73. 
S. W. Hawking, Phys. Rev. D46 (1992) 603. 

J. L. Friedmann, K. Schleich and D. N. Witt, Phys. Rev. Lett. 71 (1993) 1486. 
V. Frolov and I. D. Novikov, Phys. Rev. D48 (1993) 1607. 

M. Visser and D. Hochberg, in: Proceedings of the Haifa Workshop on the Internal 
Structure of Black Holes and Spacetime Singularities, (Haifa, Israel, 1997). 

M. S. Morris and K. S. Thorne, Am. J. Phys. 56 (1988) 395. 

M. S. Morris, K. S. Thorne and U. Yurtsever, Phys. Rev. Lett. 61 (1988) 1466. 

M. S. R. Delgaty and R. B. Mann, Int. J. Mod. Phys. D4 (1995) 231. 

P. K. F. Kuhfittig, Am. J. Phys. 67 (1999) 125. 

A. DeBenedictis and A. Das, Class. Quant. Grav. 18 (2001) 1187. 

S. Kar, N. Dadhich and M. Visser, \gT^q~c/0405103\ (2004). 

V. M. Khatsymovsky, Phys. Lett. B429 (1998) 254. 

E. Teo, Phys. Rev. D58 (1998) 024014. 

S. E. Perez Bergliaffa, K. E. Hibberd, \gr-qc/000604 1\ (2000). 
P. K. F. Kuhfittig, Phys. Rev. D67 (2003) 064015. 

M. Visser, Lorentzian Wormholes: From Einstein to Hawking. (AIP, New York, 1996). 

C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation. (Freeman, San Francisco, 
1973). 

S. Nojiri and S. D. Odintsov, Mod. Phys. Lett. A19 (2004) 627. 
G. Allemandi, A. Borowiec and M. Francaviglia, \hep-th/ 0407090 (2004). 
A. A. Starobinsky, Phys. Lett. B91 (1980) 99. 

M. C. B. Abdalla, S. Nojiri and S. D. Odintsov, Int. Conf. Math. Method. Phys. (Rio 
de Janeiro, Brazil, 2004). 

S. Carloni, P. K. S. Dunsby, S. Capozziello and A. Troisi, \gr-qc/04 lOOjo^ (2004). 
S. M. Carroll, V. Duvvuri, M. Trodden and M. S. Turner, \astro-ph/0306438\ (2003). 



11 



[26] S. Capozziello, S. Carloni and A. Troisi, Rec. Res. Dev. Astro. Astroph. RSP/AA/21 
(2003). 

X. Meng and P. Wang, Class. Quant. Grav. 21 (2004) 951. 
D. Vollick, Phys. Rev. D68 (2003) 063510. 



X. Meng and P. Wang, gr-qc/0411007, to appear in Class. Quant. Grav. (2004). 

S. Nojiri and S. D. Odintsov, Phys. Lett. B576 (2003) 5. 

U. Giinther, A. Zhuk, V. B. Bezerra and C. Romero, \hep-th/040911~2\ (2004). 

A. D. Dolgov and M. Kawasaki, Phys. Lett. B573 (2003) 1. 

T. Chiba, Phys. Lett. B575 (2003) 1. 

S. M. Carroll, A. De Felice, V. Duvvuri, D. A. Easson, M. Trodden and M. S. Turner, 
\astro-ph/041003l\ (2004). 

A. E. Dominguez and D. E. Barraco, Phys. Rev. D70 (2004) 043505. 

E. E. Flanagan, Class. Quant. Grav. 21 (2004) 3817. 
R. Dick, Gen. Rel. Grav. 36 (2004) 217. 

K. Ghoroku and T. Soma, Phys. Rev. D46 (1992) 1507. 
D. Hochberg and M. Visser, Phys. Rev. D56 (1997) 4745. 

A. DeBenedictis and A. Das, Nucl. Phys. B653 (2003) 279. 
T. A. Roman, Phys. Rev. D47 (1993) 1370. 

B. E. Taylor and W. A. Hiscock, Phys. Rev. D55 (1997) 6116. 
S. W. Kim, Grav. Cos. 6 (2000) 337. 

A. Das, A. DeBenedictis, N. Tariq, J. Math. Phys. 44 (2003) 5637. 
S. Nojiri and S. D. Odintsov, Phys. Rev. D68 (2003) 123512. 

J. P. S. Lemos, F. S. N. Lobo and S. Q. de Oliveira, Phys. Rev. D68 (2003) 064004. 

F. S. N. Lobo, Class. Quant. Grav. 21 (2004) 4811. 

F. S. N. Lobo, \gr-gc/04 10087\ (2004). 

J. L. Synge, Relativity: The General Theory (North-Holland, Amsterdam, 1964). 
W. Israel, Nuovo Cimento B44 (1966) 1. 
N. Sen, Ann. Phys. 73 (1924) 365. 

C. Lanczos, Ann. Phys. 74 518 (1924) 518. 

G. Darmois, Memorial des Sciences Mathematics XXV. (Gauthier-Villars, Paris, 1927). 



12 



